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Abstract 

Discrete maps with long-term memory are obtained from nonlinear differential equa- 
tions with Riemann-Liouville and Caputo fractional derivatives. These maps are generai- 
^n : 

izations of the well-known universal map. The memory means that their present state is 

(N ■ 

determined by all past states with special forms of weights. To obtain discrete map from 



fractional differential equations, we use the equivalence of the Cauchy-type problems and 
to the nonlinear Volterra integral equations of second kind. General forms of the universal 
maps with memory which take into account general initial conditions, for the cases of the 



^ ■ Riemann-Liouville and Caputo fractional derivatives, are suggested. 

1 Introduction 

A dynamical system consists of a set of possible states, together with a rule that determines the 
present state in terms of past states. If we require that the rule be deterministic, then we can 
define the present state uniquely from the past states. A discrete-time system without memory 
takes the current state as input and updates the situation by producing a new state as output. 
All physical classical models are described by differential or integro-differential equations, and 
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the discrete-time systems can be considered as a simplified version of these equations. A 
discrete form of the time evolution equation is called the map. Maps are important because 
they encode the behavior of deterministic systems. The assumption of determinism is that the 
output of the map can be uniquely determined from the input. In general, the present state 
is uniquely determined by all past states, and we have a discrete map with memory. Discrete 
maps are used for the study of evolution problems, possibly as a substitute of differential 
equations [U El E]- They lead to a much simpler formalism, which is particularly useful in 
simulations. The universal discrete map is one of the most widely studied maps. In this 
paper, we consider discrete maps with memory that can be used to study solutions of fractional 
differential equations [U El El E] • 

The nonlinear dynamics can be considered in terms of discrete maps. It is a very important 
step in understanding the qualitative behavior of systems described by differential equations. 
The derivatives of non-integer orders are a generalization of the ordinary differentiation of 
integer order. Fractional differentiation with respect to time is characterized by long-term 
memory effects. The discrete maps with memory are considered in [HI EJ [lOj [TTJ [12j [131 El- 
The interesting question is a connection of fractional differential equations and discrete maps 
with memory. It is important to derive discrete maps with memory from the equation of motion. 

In Ref. [14], we prove that the discrete maps with memory can be derived from differential 
equations with fractional derivatives. The fractional generalization of the universal map was 
obtained [H] from a differential equation with Riemann-Liouville fractional derivatives. The 
Riemann-Liouville derivative has some notable disadvantages such as the hyper- singular im- 
proper integral, where the order of singularity is higher than the dimension, and nonzero of the 
fractional derivative of constants, which would entail that dissipation does not vanish for a sys- 
tem in equilibrium. The desire to formulate initial value problems for mechanical systems leads 
to the use of Caputo fractional derivatives rather than Riemann-Liouville fractional derivative. 

It is possible to state that the Caputo fractional derivatives allows us to give more clear 
mechanical interpretation. At the same time, we cannot state that the Riemann-Liouville frac- 
tional derivative does not have a physical interpretation and that it shows unphysical behavior. 
Physical interpretations of the Riemann-Liouville fractional derivatives are more complicated 



than Caputo fractional derivatives. But the Riemann-Liouville fractional derivatives naturally 
appear for real physical systems in electrodynamics. We note that the dielectric susceptibil- 
ity of a wide class of dielectric materials follows, over extended frequency ranges, a fractional 
power-law frequency dependence that is called the "universal" response [T5l HE]- As was proved 
in [T71 |18], the electromagnetic fields in such dielectric media are described by differential 
equations with Riemann-Liouville fractional time derivatives. These fractional equations for 
"universal" electromagnetic waves in dielectric media are common to a wide class of materi- 
als, regardless of the type of physical structure, chemical composition, or of the nature of the 
polarizing species. Therefore, we cannot state that Riemann-Liouville fractional time deriva- 
tives do not have a physical interpretation. The physical interpretation of these derivatives in 
electrodynamics is connected with the frequency dependence of the dielectric susceptibility. As 
a result, the discrete maps with memory that are connected with differential equations with 
Riemann-Liouville fractional derivatives are very important, and these derivatives naturally 
appear for real physical systems. 

For computer simulation and physical application, it is very important to take into account 
the initial conditions for discrete maps with memory that are obtained from differential equa- 
tions with Riemann-Liouville fractional time derivatives. In Ref. [H] . these conditions are not 
obtained. In this paper, to take into account the initial condition, we use the equivalence of the 
differential equation with Riemann-Liouville and Caputo fractional derivatives and the Volterra 
integral equation. This approach is more general than the auxiliary variable method that is 
used in Ref. [H] . The proof of the result for Riemann-Liouville fractional derivatives is more 
complicated in comparison with the results for the Caputo fractional derivative. In this paper, 
we prove that the discrete maps with memory can be obtained from differential equations with 
the Caputo fractional derivative. The fractional generalization of the universal map is derived 
from a fractional differential equation with Caputo derivatives. 

The universal maps with memory are obtained by using the equivalence of the fractional 
differential equation and the Volterra integral equation. We reduce the Cauchy-type problem 
for the differential equations with the Caputo and Riemann-Liouville fractional derivatives 
to nonlinear Volterra integral equations of second kind. The equivalence of this Cauchy type 



problem for the fractional equations with the Caputo derivative and the correspondent Volterra 
integral equation was proved by Kilbas and Marzan in [191 I2Q] . We also use that the Cauchy- 
type problem for the differential equations with the Riemann-Liouville fractional derivative can 
be reduced to a Volterra integral equation. The equivalence of this Cauchy-type problem and 
the correspondent Volterra equation was proved by Kilbas, Bonilla, and Trujillo in [21, .22]. 

In Section 2, differential equations with integer derivative and universal maps without mem- 
ory are considered to fix notations and provide convenient references. In Section 3, fractional 
differential equations with the Riemann-Liouville derivative and universal maps with memory 
are discussed. In Section 4, the difference between the Caputo and Riemann-Liouville fractional 
derivatives are discussed. In Section 5, fractional differential equations with the Caputo deriva- 
tive and correspondent discrete maps with memory are considered. A fractional generalization 
of the universal map is obtained from kicked differential equations with the Caputo fractional 
derivative of order 1 < a < 2. The usual universal map is a special case of the universal map 
with memory. Finally, a short conclusion is given in section 6. 

2 Integer derivative and universal map without memory 

In this section, differential equations with derivative of integer order and the universal map 
without memory are considered to fix notations and provide convenient references. 
Let us consider the equation of motion 

oo 

Dtx{t)+KG[x{t)\Y J 5(--k)=0 (1) 

fe=i 

in which perturbation is a periodic sequence of delta-function-type pulses (kicks) following with 
period T = 2tt/u, K is an amplitude of the pulses, D\ = d 2 /dt 2 , and G[x] is some real-valued 
function. It is well-known that this differential equation can be represented in the form of the 
discrete map 

x n +i -x n = Pn+iT, p n+1 -p n = -KT G[x n ]. (2) 

Equations (|2J) are called the universal map. For details, see for example [HO [3]. 



The traditional method of derivation of the universal map equations from the differential 
equations is considered in Section 5.1 of [2J. We use another method of derivation of these 
equations to fix notations and provide convenient references. We obtain the universal map by 
using the equivalence of the differential equation and the Volterra integral equation. 

Proposition 1. The Cauchy-type problem for the differential equations 

D}x(t)=p(t), (3) 

oo 

Dlp(t) = -KG[x(t)]J2t(f-k), (4) 



fc=i 
with the initial conditions 



x(0) = x , p(0) = po (5) 

is equivalent to the universal map equations of the form 

n 

-i = x + p (n + 1)T - KT 2 J2 G W] (n + l-k), (6) 



fe=l 
n 



p n+1 =po- KT > y G[x k ] ■ (7) 

fe=i 

Proof. Consider the nonlinear differential equation of second order 

D*x{t) = G[t,x{t)], (0<t<t f ) (8) 

on a finite interval [0, tf] of the real axis, with the initial conditions 

x(0) = x , {D}x){0)=p . (9) 

The Cauchy-type problem of the form (jSJ), OH]) is equivalent to the Volterra integral equation of 

second kind 

x(t) = x +p t+ drG[r,x(r)](t-T). (10) 

Jo 

Using the function 

oo 

G[t,x(t)] = -KG[x(t)]J2^(^-k), 



fe=i 
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for nT <t < (n+ 1)T, we obtain 

n 

x(t) = x + Po t - KT^G[x{kT)} (t - kT). (11) 

fe=i 

For the momentum p(t) = Dixit), equation ( 1TT1) gives 

n 

P (t)=p -KTj2G[x(kT)]. (12) 

fc=i 

The solution of the left side of the (n+ l)th kick 

x n+1 = x(t n+1 - 0) = lim x(T(n + 1) - e), (13) 

Pn+i = p(t n+ i - 0) = lim p(T(n + 1) - e), (14) 

where t n+ i = (n + 1)T, has the form (jBj) and (j7|). 
This ends the proof. D 

Remark 1. 

We note that equations (jSJ) and (|7|) can be rewritten in the form (J2J). Using equations © and 
(|II), the differences x n+ i — x n and p n +i — Pn give equations (J2J) of the universal map. 

Remark 2. 
If G[x] = —x, then equations ([2]) give the Anosov-type system 

X n+ i ~X n = Pn+lT, Pn+1 ~Pn = KTx n . (15) 

For G[x] = sin(x), equations (J2J) are 

x n+1 -x n = Pn+iT, p n+1 -p n = -KT sin(x n ). (16) 

This map is known as the standard or Chirikov- Taylor map PQ. 



3 Riemann-Liouville fractional derivative and universal 
map with memory 

In this section, we discuss nonlinear differential equations with the left-sided Riemann-Liouville 
fractional derivative $Df defined for a > by 

1 d n /"* rlYWr 

D?x(t) = D n t I^ a x(t) = TT7 c-j- / 1 \i , (n-l<a<n), (17) 

where D™ = d n /dt n , and qI" is a fractional integration [U [7J E]. 
We consider the fractional differential equation 

D?x(t) =G[t,x(t)], (18) 

where G[t,x(t)} is a real-valued function, < n — 1 < a < n, and t > 0, with the intial 
conditions 

( D?- k x)(0+) = a*, fc = l,...,n. (19) 

The notation ( -D" _fc a;)(0-|-) means that the limit is taken at almost all points of the right-sided 
neighborhood (0, + e), e > 0, of zero as follows 

(oA Q ~^)(0+) = t hm + D?- k x(t), (k = 1, ...,n- 1), 
( D?~ n x)(0+) = lim I?- a x(t). 



The Cauchy-type problem (1181) and ([191 can be reduced to the nonlinear Volterra integral 
equation of second kind 

where £ > 0. The result was obtained by Kilbas, Bonilla, and Trujillo in J2TJ[22]. For a = n = 2, 
equation (1201) gives (TTUj) . 

The Cauchy-type problem ffTB"]) and ffT9|) and the nonlinear Volterra integral equation ( 120]) 
are equivalent in the sense that, if x{t) e L(0,t/) satisfies one of these relations, then it also 



satisfies the other. In [2TJ [22] (see also Theorem 3.1. in Section 3.2.1 of [7]), this result is 
proved by assuming that the function G[t, x] belongs to L(0, tf) for any x G W C R. 

Let us give the basic theorem regarding the nonlinear differential equation involving the 
Riemann-Liouville fractional derivative. 

Kilbas-Bonilla-Trujillo Theorem. Let W be an open set in R and let G[t,x], where 
t G (0, tf] and x G W , be a real-valued function such that G\t,x] G L(0,tf) for any x G W . Let 
x(t) be a Lebesgue measurable function on (0,tf). If x(t) G L(0,tf), then x(t) satisfies almost 
everywhere equation l[T8\) and conditions $TQ) if, and only if, x(t) satisfies almost everywhere 
the integral equation (05J). 



Proof. This theorem is proved in J2~T||2"2"] (see also Theorem 3.1. in Section 3.2.1 of [7]). D 

In Ref. [H] we consider a fractional generalization of equation (0Q) of the form 

D?x(t) + KG[x(t)]J2^(f -k)=0, (1< a < 2), (21) 

fc=i 

where t > 0, and $Df is the Riemann-Liouville fractional derivative defined by (fTTj) . Let us 
give the following theorem for equation (!2~T|) . 

Proposition 2. The Cauchy-type problem for the fractional differential equation of the 
form ( [HP with the initial conditions 

( Dr 1 x)(0+) = Cl , ( Dr 2 x)(0+) = ( / t 2 - a x)(0+) = c 2 (22) 

is equivalent to the equation 

r(a) r(a-l) r(a) ^ 

where nT < t < (n + 1)T. 



Proof. Using the function 



oo 



G[t,z{t)] = -KG[z]Y l S(f- k )> ( 24 ) 

fe=i 

equation 1 J2"T|) has the form of ( fT8|) with the Riemann-Liouville fractional derivative of order a, 

where 1 < a < 2. It allows us to use the Kilbas-Bonilla-Trujillo theorem. As a result, equation 

(|2ip with initial conditions f lT9|) of the form (122]) is equivalent to the nonlinear Volterra integral 

equation 

is °° rt 



J2 / drG^rJKt-r)- 1 *^-*), (25) 

fc=i ^° 



r(a) r(a-l) r(a) 

where t > 0. If nT < £ < (n + 1)T, then the integration in (1251) with respect to r gives ( [231 . 
This ends the proof. □ 

To obtain equations of discrete map a momentum must be defined. There are two possiblities 
of defining the momentum: 

p(t) = D?- l x(t), P(t) = D]x{t). (26) 

Let us use the first definition. Then the momentum is defined by the fractional derivative 
of order a — 1. Using the definition of the Riemann-Liouville fractional derivative ( !T7|) in the 
form 

D?x(t) = D 2 Q I 2 t ~ a x{t), (1< a < 2), (27) 

we define the momentum 

pit) = oDf^xCt) = — ^ r4 / , ^^ 1 , (K « < 2), (28) 

Fy ' l w r(2 - a) dt J (t-r)"" 1 v ~ ; y J 

where x(t) is defined for r G (0,t). Then 

D?x(t) = D]p{t), (Ka<2). (29) 

Using momentum p[t) and coordinate x(t), equation ( 12T1) can be represented in the Hamiltonian 
form 

D?- 1 x{t)=p{t) 1 (30) 



x t 



D]p{t) = -KG[x{t)}^{- f -k). (31) 

fc=i 

Proposition 3. The Cauchy type problem for the fractional differential equations of the 
form /[3ty) and l[3l\) with the initial conditions 

( Dr 1 x)(0+) = Cl , ( D?- 2 x)(0+) = ( / t 2 - Q x)(0+) = c 2 (32) 

is equivalent to the discrete map equations 

„ rpa-l rp a -2 KT a n 

*n + i = %r^(n + ir 1 + ^ V An + l) a ~ 2 -YT^H G N (" + X - fc ) a ~ 1 ' ( 33 ) 

r(a) r(a-l) r(a) ^ 

n 

Pn+1 = Cl -irT^G[x fc ]. (34) 

fe=i 

Proof. We use Proposition 2 to prove this statement. To obtain an equation for the 
momentum (1281) . we use the following fractional derivatives of power functions (see Section 2.1 

inH), 

a Df{t-af- 1 = J} P) (t-af- 1 -*, a>0, (3 > 0, t > a, (35) 

l [p — a) 

D?t a ~ k = 0, k = l,...,n, n-l<a<n. (36) 



O^t 



These equations give 



and 



Dff~ l = r(a) 



Dft a - 2 = 0, 



a A Q (t -<=!». 



We note that equation (I23p for x(r) can be used only if r £ {nT, t), where nT < t < (n+l)T. 
The function x(r) in the fractional derivative qD^ of the form fl28l) must be defined for all 
r £ (0, t). We cannot take the derivative Df of the functions (r — fcT) a_1 that are defined 
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for r G (kT,t). In order to use equation (123]) on the interval (0, £), we must modify the sum in 
equation ( 123]) by using the Heaviside step function. Then equation (123]) has the form 

T(a) T(a- 1) r(a) ^ 

where r G (0, £). Using the relation 

Q D?((t - a) a ~ l e(t - a)) = a D?(t - a)*" 1 = I», (38) 

equations (128]) and ( 137]) give 

p(t) =ci-KTJ2 G[x(kT)}, (39) 

fe=i 

where nT < t < (n + 1)T. Then the solution of the left side of the (n + l)-th kick 

p n+1 = Cl - KT J2 G[x k ]. (40) 

fe=i 

As a result, we obtain a universal map with memory in the form of equations ( 133]) and (j3i 
This ends the proof. □ 

Remark 3. 

For a = n = 2 equations ( 133]) and ( 134]) give the usual universal map ([6]) and (J7J. 

Remark 4. 
We note that the map (133]) and ( 1341 with 

ci = pi, c 2 = 

was obtained in [T3] in the form 

T 



-.Q-l ™ 



^n+1 



v y fc=i 



5^p fc+ iV;(n-A; + l), (41) 



Pn+1 = Pn -KTG(x n ), (Ka<2), (42) 
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where p\ — C\, and the function V a (z) is defined by 

V a (z) = z a ' 1 -(z-ir-\ (z>l). (43) 

In Ref. [H], we obtain these map equations by using an auxiliary variable £(£) such that 

The nonlinear Volterra integral equations and the general initial conditions ( 13 2 p are not used 
in [H]. In the general case, the fractional differential equation of the kicked system (121]) is 
equivalent to the discrete map equations 

Xn+1 = Y(aj £ P k + lV ^ n ~ k + l ) + rll-l) ^ + ir " 2 ' (44) 

Pn+1 = Pn -KTG(x n ), (Ka<2), (45) 

where pi = c\. Here we take into account the initial conditions (!32|l . The second term of the 
right-hand side of equation (|4"4"1) is not considered in [Hj. Using — l<a; — 2<0, we have 

lim (n + l) a " 2 = 0. 

n— >oo 

Therefore, the case of large values of n is equivalent to c<i = 0. 

Let us give the proposition regarding the second definition of the momentum p(t) = D]x(i). 

Proposition 4. The Cauchy-type problem for the fractional differential equations 

D]x(t)=p(t), (46) 

oo 

D?x{t) = -KG[x(t)}J2s( f -k), (1< a < 2), (47) 



wt/i i/ie initial conditions 

( o Dr 1 x)(0+) = c l7 ( 0J Dr 2 x)(0+) = ( / 2 - a x)(0+) = c 2 (4? 
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is equivalent to the discrete map equations 



:r 



n+l 



(n + l)"" 1 + J—T-dn + ir 2 " w-y E G N(™ + 1 - fc) 0-1 - (49) 



Pn+l 



I» " ' Y{a-\y ' I» ^ 

c i T "" 2 /.. , i**-a , c 2 (« - 2)T°-*, , __ 3 ffT- 1 



( w + 1 )°" 2 + rfa-1) ( " +ir " 3 ~lY^iy S GN(nH-l-fc)- 2 . (50) 



i>-i) v y r(a-i) A ..., 



Proof. We define the momentum 

p(t) = D]x{t). 
If nT < t < (n + 1)T, then the differentiation of (123]) with respect to £ gives 

Here we use the relation 

T(a) = (a-l)T(a-l), (1 < a < 2). 

Using equations ([23]) and (15 ip . we can obtain the solution of the left side of the (n + l)-th kick 
(fT3"|) and (|T4|) . As a result, we have equations (|4"9~]l and (1501 . 
This ends the proof. □ 

Remark 5. 

Equations ( f49l) and ( [501) describe a generalization of equations (jEJ) and (CO). If a = n = 2, and 
c 2 = x , Ci = p 0) then equation (1491) gives dHj) and (J7J). 

Remark 6. 
In equations ( 1501) and ( J5T1) . we can use 

c 2 (a - 2) _ c 2 

r(a- 1) ~ r(a-2) 

for 1 < a < 2. 
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Remark 7. 

If we use the definition pit) = D\x{t), then the Hamiltonian form of the equations of motion 
will be more complicated than ( 130]) and ( 13TJ) since 

D\ lt a x(t) ± ,P t - a D\x{t). 

Remark 8. 

Note that we use the usual momentum pit) = D\x{t). In this case, the values c\ and C2 are not 
connected with p(Q) and x(0). If we use the momentum pit) = qD^~ x(t), then c\ = p(0). 

4 Riemann-Liouville and Caputo fractional derivatives 

In Ref. [H] we consider nonlinear differential equations with Riemann-Liouville fractional 
derivatives. The discrete maps with memory are obtained from these equations. The problems 
with initial conditions for the Riemann-Liouville fractional derivative are not discussed. 

The Riemann-Liouville fractional derivative has some notable disadvantages in applications 
in mechanics such as the hyper- singular improper integral, where the order of singularity is 
higher than the dimension, and nonzero of the fractional derivative of constants, which would 
entail that dissipation does not vanish for a system in equilibrium. The desire to use the usual 
initial value problems for mechanical systems leads to the use of Caputo fractional derivatives 
[7J E] rather than the Riemann-Liouville fractional derivatives. 

The left-sided Caputo fractional derivative [23j I2U [25j [7J of order a > is defined by 



oVtI[t) ~T(n-a)J (t-r)"-'* 1 



D tf(t) = ^7— — r / n :C> = oir a Dtf(t), (52) 



where n — 1 < a < n, and qI" is the left-sided Riemann-Liouville fractional integral of order 
a > that is defined by 

This definition is, of course, more restrictive than the Riemann-Liouville fractional derivative 
[U [7] in that it requires the absolute integrability of the derivative of order n. The Caputo 
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fractional derivative first computes an ordinary derivative followed by a fractional integral to 
achieve the desire order of fractional derivative. The Riemann-Liouville fractional derivative is 
computed in the reverse order. Integration by part of (152]) will lead to 

?£>?*(*) = D?x(t) - J2 v(1 f a , u xW(0). (54) 

It is observed that the second term in equation f l54l) regularizes the Caputo fractional derivative 
to avoid the potentially divergence from singular integration at t = 0. In addition, the Caputo 
fractional differentiation of a constant results in zero 

o c A Q c = o. 

Note that the Riemann-Liouville fractional derivative of a constant need not be zero, and we 
have 



D?C = — -C. 



r(T^a) 

If the Caputo fractional derivative is used instead of the Riemann-Liouville fractional deriva- 
tive, then the initial conditions for fractional dynamical systems are the same as those for the 
usual dynamical systems. The Caputo formulation of fractional calculus can be more applicable 
in mechanics than the Riemann-Liouville formulation. 

5 Caputo fractional derivative and universal map with 

memory 

In this section, we study a generalization of differential equation ([1]) by the Caputo fractional 
derivative. The universal map with memory is derived from this fractional equation. 
We consider the nonlinear differential equation of order a, where <n — 1 < a < n, 

oD a x(t) = G[t, x{t)\, (0 < t < t f ), (55) 

involving the Caputo fractional derivative $Df on a finite interval [0, tf] of the real axis, with 
the initial conditions 

(D*x)(0) =c fc , k = 0,...,n-l. (56) 

15 



Kilbas and Marzan (T9J EH] proved the equivalence of the Cauchy-type problem of the form 
fl55|) . fl56l) and the Volterra integral equation of second kind 

n— 1 _. »t 

x ^ = J2 c i tk + f^) dT G[r ' x(r)] (t - r)Q_1 (57) 

fc=0 ' I J "'O 

in the space C n_1 [0,£/]. For a = n = 2 equation f )57|) gives ffTUj) . 

Let us give the basic theorem regarding the nonlinear differential equation involving the 
Caputo fractional derivative. 



Kilbas-Marzan Theorem. The Cauchy-type problem / T55]) and l[5b)) and the nonlinear 
Volterra integral equation f57| ) are equivalent in the sense that, if x(t) 6 C[0,t/] satisfies one 
of these relations, then it also satisfies the other. 

Proof. In [T9l [20] (see also [7], Theorem 3.24.) this theorem is proved by assuming that a 
function G[t, x] for any x G W C M. belong to C 7 (0,t/) with < 7 < 1, 7 < a. Here C 7 (0,t/) 
is the weighted space of functions f[t] given on (0,tf], such that t 7 /[t] £ C(0,tf). This ends 
the proof. □ 



We consider the fractional differential equation of the form 

fe=i 
where $Df is the Caputo fractional derivative, with the initial conditions 



00 
%D?x{t) +KG[x(t)]J2t(f ~ k) = 0, (1< a < 2), (5J 



x(0) = x , ( J D 1 a;)(0)=po- (59) 

Using p(t) = D\x{t), equation (|58|) can be rewritten in the Hamilton form. 

Proposition 5. TTie Cauchy-type problem for the fractional differential equations 

Dlx(t)=p(t), (60) 
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x t 



^ryt) = -KG[x(t)]^5(--fc), (Ka<2), (61) 

fe=i 

mi/i t/ie initial conditions 

x(0) = x , p(O) = po (62) 

zs equivalent to the discrete map equations 

x n+ i =x + p (n + 1)T - — - V (ti + 1 - k) a ' x G[x k i (63) 

Pn+1 =Po ~ T(a-l) T,( n + 1 - k) a - 2 G[x k }. (64) 

^ ' fc=i 

Proof. We use the Kilbas-Marzan theorem with the function 

fe=i 
The Cauchy-type problem (1551) and (159]) is equivalent to the Volterra integral equation of second 



oo 

G[t,x(t)] = -^G[x(t)]^5(--fc). 



kind 

T^ °° pt _L 

x (t) = x +p t-—-J2 dT(t~rr- l G[x(r)}5^ f -ky (65) 

in the space of continuously differentiable functions x(t) 6 C 1 [0,t/]. 
If nT < t < (n + 1)T, then equation (1631) gives 

a#) = Xo + Po t _ V (t - k:T) a - 1 G[x(kT)). (66) 

We define the momenta 

p{t) = D]x{t). (67) 

Then equations (166]) and (IB"?]) give 

P(*) = Pb - tv TT E (* - kT) a - 2 G[x(kT)}, (nT<t<(n + 1)T), (68) 

where we use iYa) = (a — l)r(a — 1). 
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The solution of the left side of the (n + l)-th kick (JT3ll and (1141) can be represented by 
equations fl63|) and f )64l) . where we use the condition of continuity x(t n + 0) = x(t n — 0). 
This ends the proof. □ 

Remark 9. 

Equations ( )63|) and ( )64|) define a generalization of the universal map. This map is derived from 
a fractional differential equation with Caputo derivatives without any approximations. The 
main property of the suggested map is a long-term memory that means that their present state 
depends on all past states with a power-law form of weights. 

Remark 10. 
If a = 2, then equations f|63|) and (|64"|) give the universal map of the form ([6]) and ((7|) that 
is equivalent to equations (|2]). As a result, the usual universal map is a special case of this 
universal map with memory. 

Remark 11. 
By analogy with Proposition 5, it is easy to obtain the universal map with memory from 
fractional equation ( 1581) with a > 2. 



6 Conclusion 

The suggested discrete maps with memory are generalizations of the universal map. These maps 
describe fractional dynamics of complex physical systems. The suggested universal maps with 
memory are equivalent to the correspondent fractional kicked differential equations. We obtain 
a discrete map from fractional differential equation by using the equivalence of the Cauchy-type 
problem and the nonlinear Volterra integral equation of second kind. An approximation for 
fractional derivatives of these equations is not used. 

It is important to obtain and to study discrete maps which correspond to the real physical 
systems described by the fractional differential equations. Media with memory in mechanics and 
electrodynamics, we can consider viscoelastic and dielectric materials as a media with memory. 
We note that the dielectric susceptibility of a wide class of dielectric materials follows, over 
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extended frequency ranges, a fractional power-law frequency dependence that is called the 
"universal" response [15} IT6] . As was proved in [TTt IT8] . the electromagnetic fields in such 
dielectric media are described by differential equations with fractional time derivatives. These 
fractional equations for electromagnetic waves in dielectric media are common to a wide class of 
materials, regardless of the type of physical structure, chemical composition, or of the nature of 
the polarizing species, whether dipoles, electrons or ions. We hope that it is possible to obtain 
the discrere maps with memory which correspond to the real dielectric media described by the 
fractional differential equations. 

Acknowledgments 

This work was supported by the Office of Naval Research, Grant No. N00014-02-1-0056 and 
Rosnauka No.02.740.11.0244. 

References 

[1] B.V. Chirikov, "A universal instability of many dimensional oscillator systems" Physics 
Reports 52(5), (1979) 263-379. 

[2] R.Z. Sagdeev, D.A. Usikov, G.M. Zaslavsky, Nonlinear Physics: From the Pendulum to 
Turbulence and Chaos (Harwood Academic Publishers, New- York, 1988). 

[3] G.M. Zaslavsky, Hamiltonian Chaos and Fractional Dynamics (Oxford University Press, 
Oxford, 2005). 

[4] S.G. Samko, A. A. Kilbas, O.I. Marichev, Fractional Integrals and Derivatives Theory and 
Applications (Gordon and Breach, New York, 1993). 

[5] K. Miller, B. Ross, An Introduction to the Fractional Calculus and Fractional Differential 
Equations (Wiley, New York, 1993). 

[6] I. Podlubny, Fractional Differential Equations (Academic Press, San Diego, 1999). 

[7] A. A. Kilbas, H.M. Srivastava, J.J. Trujillo, Theory and Application of Fractional Differ- 
ential Equations (Elsevier, Amsterdam, 2006). 

19 



[8] A. Fulinski, A.S. Kleczkowski, "Nonlinear maps with memory" Physica Scripta 35(2), 

(1987) 119-122. 
[9] E. Fick, M. Fick, G. Hausmann, "Logistic equation with memory" Physical Review A 

44(4), (1991) 2469-2473. 
[10] K. Hartwich, E. Fick, "Hopf bifurcations in the logistic map with oscillating memory" 

Physics Letters A 177(4/5), (1993) 305-310. 
[11] M. Giona, "Dynamics and relaxation properties of complex systems with memory" Non- 
linearity 4(3), (1991) 991-925. 
[12] J.A.C. Gallas, "Simulating memory effects with discrete dynamical systems" Physica A 

195(3/4), (1993) 417-430; "Erratum" Physica A 198(1/2), (1993) 339-339. 
[13] A. A. Stanislavsky, "Long-term memory contribution as applied to the motion of discrete 

dynamical system" Chaos 16(4), (2006) 043105. 
[14] V.E. Tarasov, G.M. Zaslavsky, "Fractional equations of kicked systems and discrete maps" 

Journal of Physics A 41(43), (2008) 435101. ( arXiv: 1107.39531 

[15] A.K. Jonscher, Universal Relaxation Law, (Chelsea Dielectrics Press, London, 1996). 

[16] A.K. Jonscher, "Dielectric relaxation in solids" J. Physics D Appl. Phys. 32 (1999) R57- 
R70. 

[17] V.E. Tarasov, "Universal electromagnetic waves in dielectric" Journal of Physics: Con- 
densed Matter. Vol.20. No. 17. (2008) 175223. ( arXiv:0907.2l63JI 

[18] V.E. Tarasov, "Fractional equations of Curie- von Schweidler and Gauss laws" Journal of 
Physics: Condensed Matter. Vol.20. No. 14. (2008) 145212. (larXiv:0907.1837D 



[19] A. A. Kilbas, S.A. Marzan, "The Cauchy problem for differential equations with fractional 
Caputo derivative" Doklady Mathematics 70(3), (2004) 841-845; Translated from Doklady 
Akademii Nauk 399(1), (2004) 7-11. (in Russian). 

[20] A. A. Kilbas, S.A. Marzan, "Nonlinear differential equations with the Caputo fractional 
derivative in the space of continuously differentiable functions" Differential Equations 
41(1), (2005) 84-89; Translated from Different sialnye Uravneniya 41(1), (2005) 82-86. 
(in Russian). 

20 



[21] A. A. Kilbas, B. Bonilla, J.J. Trujillo, "Nonlinear differential equations of fractional order 
is space of integrable functions" Doklady Mathematics 62(2), (2000) 222-226; Translated 
from Doklady Akademii Nauk 374(4), (2000) 445-449. (in Russian). 

[22] A. A. Kilbas, B. Bonilla, J.J. Trujillo, "Existence and uniqueness theorems for nonlinear 
fractional differential equations" Demonstratio Mathematica 33(3), (2000) 583-602. 

[23] M. Caputo, "Linear models of dissipation whose Q is almost frequency independent. Part 
II" Geophysical Journal of the Royal Astronomical Society 13(5), (1967) 529-539. 

[24] M. Caputo, F. Mainardi, "A new dissipation model based on memory mechanism" Pure 
and Applied Geophysics 91(1), (1971) 134-147. 

[25] A. Carpinteri, F. Mainardi (Eds), Fractals and Fractional Calculus in Continuum Mechan- 
ics (Springer, Wien, 1997). 

[26] V.E. Tarasov, Fractional Dynamics: Application of Fractional Calculus to Dynamics of 
Particles, Fields and Media (Springer, HEP, 2011) 504 pages. 



21 



